Tunable dipolar resonances and Einstein-de Haas effect in a 87 Rb atoms condensate 
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We study a spinor condensate of 87 Rb atoms in F = 1 hyperfine state confined in an optical dipole 
trap. Putting initially all atoms in irif = 1 component we observe a significant transfer of atoms to 
other, initially empty Zeeman states exclusively due to dipolar forces. Because of conservation of a 
total angular momentum the atoms going to other Zeeman components acquire an orbital angular 
momentum and circulate around the center of the trap. This is a realization of Einstein-de Haas 
effect in a system of cold gases. We show that the transfer of atoms via dipolar interaction is possible 
only when the energies of the initial and the final sates are equal. This condition can be fulfilled 
utilizing a resonant external magnetic field, which tunes energies of involved states via the linear 
Zeeman effect. We found that there are many final states of different spatial density which can 
be tuned selectively to the initial state. We show a simple model explaining high selectivity and 
controllability of weak dipolar interactions in the condensate of 87 Rb atoms. 



I. INTRODUCTION 

Experimental achievement of Bose-Einstein condensa- 
tion in a gas of 52 Cr atoms [l|, [|[ has launched a huge 
interest in studying properties of ultracold dipolar sys- 
tems [!, |3] • Although polar molecules with their large 
electric dipole moments seem to be perfect candidates to 
investigate dipolar effects such systems have not been 
condensed so far. Also the necessity of an external 
field inducing the dipole in these systems might suppress 
some intrinsic interaction effects (for instance, the elec- 
tric dipolc-dipole interactions do not conserve the total 
angular momentum, i.e., the sum of spin and orbital an- 
gular momentum). On the other hand, experiments with 
chromium condensates already showed spectacular fea- 
tures characteristic for dipolar interactions. For exam- 
ple, in Ref. [|[ the suppression (and even complete in- 
hibition) of inversion of ellipticity during expansion of a 
cloud of chromium atoms was observed. In Ref. @, on 
the other hand, the collapse dynamics of 52 Cr conden- 
sate was studied. The scattering length was decreased by 
means of a Feshbach resonance below the critical value 
above which the condensate is stable. Then the system 
was hold in a trap for some time and eventually the trap 
was switched off. The specific cloverleaf patterns for the 
density caused by anisotropic dipole-dipole interactions 
were observed when the system was imaged after time of 
flight. 

Some authors suggested, however, that the magnetic 
dipolar interactions might already lead to observable ef- 
fects in condensates of alkali atoms j7j-ll3|. whose mag- 
netic dipole moment is an order of magnitude lower than 
that of chromium atoms. Indeed, in experimental work 
of Ref. [lj] it is demonstrated that a spin-1 87 Rb spinor 
condensate exhibits dipolar properties. The authors of 



[Tij observe the spontaneous decay of helical spin tex- 
tures toward a spatially modulated structure of spin do- 
mains and relate this effect to the presence of magnetic 
dipolar interactions in the system. This claim is sup- 
ported by an observation that the modulated phase is 
suppressed when the dipolar interactions are reduced. 
There are, however, some difficulties in understanding 
the results of this experiment within the frame of mean- 
field and Bogoliubov theories at zero temperature fTo] ]. 

One of the most spectacular manifestations of dipo- 
lar interactions is the Einstein-de Haas effect [l|| stud- 
ied recently theoretically in gaseous systems of chromium 
(l7l [l8j and rubidium condensates 0] . Interestingly, this 
phenomenon might occur even in very weak dipolar sys- 
tems. To observe the Einstein-de Haas effect in alkali 
atoms condensates one has to tune an external magnetic 
field to some resonant value or for a fixed magnetic field 
to modify appropriately trapping parameters. Only when 
the Zeeman energy fits the energy gap between different 
single-particle spatial states associated with mj? = 1 and 
mp = spin states, the transfer of atoms is resonantly 
amplified. 

We note that the resonant magnetic field is typically 
of the order of tens or hundreds micro-Gauss making 
the observation of the Einstein-de Haas effect difficult at 
present. Since dipole-dipole interactions are very weak 
the resonant curves are also very narrow. This means 
that experimental realization needs high precision. On 
the other hand, it guarantees that resonances are highly 
selective. By choosing an appropriate value of the mag- 
netic field one can tune the transition of atoms to par- 
ticular spatial state. Indeed, controlling dipolar interac- 
tions is the crucial point in working with dipolar systems. 
Such a control has been recently imposed in chromium 
condensate 113. It was shown that the external static 
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magnetic field strongly influences the dipolar relaxation 
rate - there exists a range of magnetic field intensities 
where this relaxation rate is strongly reduced allowing 
for the accurate determination of S = 6 scattering length 
for chromium atom. In Rcf . [20] , on the other hand, two- 
dimensional optical lattice and a static magnetic field are 
used to control the dipolar relaxation into higher lattice 
bands. In this work an evidence for the existence of the 
relaxation threshold with respect to the intensity of the 
magnetic field is given. As the authors of Ref. [2(| claim 
such an experimental setup might lead to the observation 
of the Einstein-de Haas effect. 

This paper is organized as follows. In Sec. [TT] we de- 
scribe the system under consideration. In particular, 
we discuss the contact and dipolar interactions (we re- 
fer the reader to Appendix [Al for further details). Sec. 
IIIII presents numerical results regarding the resonances 
and the Einstcin-de Haas effect in a rubidium conden- 
sate whereas Sec. IIVI introduces a simple model explain- 
ing the origin of the resonances as well as properties of 
spatial states populated via dipolar interactions. We end 
with conclusions in Sec. [V] Moreover, the Appendix [B] 
offers a lot of technical details implemented in numeri- 
cal procedure used to obtain the resultsjsresented in this 
paper as well as in our previous work 



II. DESCRIPTION OF THE SYSTEM 

We investigate a spinor condensate of atoms in the 
F = 1 hyperfine state. In addition to the dominating bi- 
nary contact interactions, we consider long-range dipolar 
magnetic interactions. In the formalism of second quan- 
tization, the Hamiltonian of the system is given by 
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with repeated indices being summed over the values +1, 
0, and —1. The field operator V>i(r) (V'J(t')) annihi- 
lates (creates) an atom in the hyperfine state \F = 
at point r. The first row in (pj represents the single- 
particle Hamiltonian that consists of the kinetic energy 
part and the trapping potential assumed to be indepen- 
dent of the hyperfine state. The first term in the second 
row describes the interaction with the magnetic field B 
which is in our case always directed along the z-axis. 
Magnetic moment of the atom fj, if proportional to the 
dimensionless spin-I matrices F, fj, = fj,F. Effective 
magnetic moment of the atom fj, differs from the Bohr 
magneton by Lande giromagnetic factor 7, fi = 7/Xb. 



The terms with coefficients cq and ci describe the spin- 
independent and spin-dependent parts of the contact in- 
teractions, respectively. Coefficients cq and C2 can be 
expressed with the help of the scattering lengths eto and 
02 which determine the collision of atoms in a channel of 
total spin and 2. One has cq = 4-7rft 2 (ao + 202) /3m and 
C2 = Airh 2 (a2 — ao)/3m (m is the mass of an atom) [2lj . 
where ao = 5.387 nm and 0,2 = 5.313 nm (22j. Finally, 
the last term represents the magnetic dipolar interac- 
tions and originates from the interaction energy of two 
magnetic dipole moments [i 1 and fj, 2 
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where r and r' denote the positions of the dipoles. 
The equation of motion reads 



= {% 





(3) 



where % sp = — + V ex t — fimpB represents the 

single-particle part (including the interaction with the 
external magnetic field) whereas H c c and H d originate 
from the two-particle interactions and are discussed in 
detail in the Appendix [SJ Here we just state that the 
diagonal part of T~L C is responsible for elastic collisions 
whereas other elements of H c allow for the change of 
spin projections of individual atoms preserving, however, 
the projection of total spin. 

The dipolar part of Eq. ([3]), H d , is more complex. To 
better understand the dipolar processes it is convenient 
to rewrite the dipolar interactions in the following form 

m 
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where Y2 fl (f) (with r denoting a unit vector in the di- 
rection of relative position of two atoms) is a spherical 
harmonics of rank-2 and £2^ defined as 



£2,0= -\ MzF 2z -F 1 -F 2 /3) 



S 2 ,±i = ±-(F u F 2 ± + F 1± F 2z ) 
£2, ±2 = — —Fr±F%± 



(5) 



is a rank-2 spherical tensor built of atomic spin operators 
(including raising and lowering operators -Fi± and F%±). 
It is clear from (|4]) that when two atoms interact the 
total spin projection Mp (as well as the total spin itself) 
can change at most by 2 whereas the spin projection of 
individual atoms (see (O) changes maximally by 1 not 
by 2. Therefore, the atom can not be transferred directly 
from the mp = 1 torn^ = —1 component (populating 
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the mp = — 1 state is a second order process). When 
AMp = ±2, the last row in ((5|) shows that the projection 
of the spin of each atom changes by +h or —h, i.e., both 
atoms initially in the same state go simultaneously to 
the nearest (in a sense of magnetic number mp) state or 
when atoms arc in different but neighboring components 
they are transferred to the states shifted in number mi? 
by +1 or —1. In addition to the processes just described 
there are the atomic collisions in which only one atom is 
transferred to the other Zccman state (AMp = ±1, see 
the middle row in ([5|)) or the collisions that do not change 
the spin projection (AMp = 0, the first row in $5§). In 
the latter case, however, still the process transferring one 
atom from m^ = 0to?7iF = l and the other from mp = 
to mp = — 1 states is allowed just like it is allowed 
for contact interactions. According to ([J) in each case 
the change of the total spin projection is accompanied 
by an appropriate change of a relative orbital angular 
momentum of colliding atoms. 

Hence, the dipolar interactions do not conserve the 
projection of total spin of two interacting atoms. Nei- 
ther the projection of total orbital angular momentum 
is preserved (see ©)• However, the dipolar interactions 
couple the spin and the orbital motion of atoms (with 
regard to any axis) as revealed by the last relation in © 

[V d , F lz + F 2z ] 0, {V d , L u + L 2z ] ? 0, 
[V d ,L lz + L 2z + F lz + F 2z ] =0. (6) 

Therefore, going to mp = and mp = — 1 states atoms 
acquire the orbital angular momentum and start to cir- 
culate around the center of the trap. This is the essence 
of the famous Einstein-de Haas effect and the discussion 
in Sec. IIIII suggests its possible realization in cold gases. 

III. DIPOLAR RESONANCES AND THE 
EINSTEIN-DE HAAS EFFECT 

In this Section we investigate the properties of dipolar 
spinor condensate of 87 Rb atoms in the F = 1 hypcrfinc 
state (7 = 1/2) by solving the mean- field version of Eq. 


lh i l^|= (n sp + n c + n d ) ( to ) (7) 

with ipi being the wave function (referred sometimes also 
to as an order parameter) of the i-th spin component. 
Usually, dipolar properties of rubidium condensates are 
neglected. This is because the magnetic moment of ru- 
bidium atoms is small. Comparing typical energy re- 
lated to the dipolar interactions (given by fJ?n, where 
/j is the magnetic moment of 87 Rb atom and n is the 
atomic density) to the characteristic contact interactions 
energy gn (g determines the strength of the contact inter- 
actions) one obtains for the maximally stretched F = 1 
hypcrfinc state (i.e., for the spin-polarized case where 



g = 47rft 2 a 2 /m) of 87 Rb the ratio of the order of 10~ 4 . 
This is two orders of magnitude lower than the corre- 
sponding ratio calculated for condensate of chromium 
atoms which is commonly considered as a dipolar con- 
densate. In this Section we show, however, that for a 
rubidium condensate there exist resonances which sig- 
nificantly enhance the effect of dipolar interactions and 
might lead to the manifestation of various properties in 
as strong way as in the case of chromium condensate. 

The numerical experiment, we perform, consists in 
preparing the condensate initially m my = I Zeeman 
state (all atomic magnetic moments aligned along the 
magnetic field). It is done by evolving the mean-field 
version of Eq. <j3j) in imaginary time at the presence of 
external magnetic which is large enough to enforce al- 
most all atoms to populate mp = 1 component. Next, 
we change both the direction and the value of the mag- 
netic field. What we usually observe is the transfer of 
small number of atoms to mp = and rrif = -1 com- 
ponents. Wc would like to stress that in our numerical 
experiment the initial condition excludes the short-range 
spin dynamics and the depletion of mf = 1 state occurs 
only due to dipolar interactions. 

However, there exist the particular values of the mag- 
netic field (for example, B w — O.llmG for the axi- 
ally symmetric trap with u z = 2n x 100 Hz and lo t = 
2ir x 400 Hz and initial number of atoms TV = 5 x 10 4 , 
as it is shown in Fig. [1]) when the observed transfer is 
large - of the order of 10% of initial number of atoms in 
mp = 1 state. The maximal transfer happens after the 
time of the order of 100 ms (sec Fig. \2§ and is consis- 
tent with the characteristic time scale related to dipolar 
interaction, which is h//i 2 n. This formula tells us that 
increasing the atomic magnetic moment (for example, by 
going from rubidium to chromium condensate) the time 
needed for reaching the maximal transfer gets shorter. 
Indeed, for chromium condensate (with 12 times larger 
the magnetic moment in comparison with 87 Rb atoms in 
F = 1 state) this time is of the order of 1 ms [I?} . Similar 
behavior is visible for larger magnetic fields. 

An interesting observation is made when looking at 
the phase and the density of spinor components for all 
resonances we found (Figs. [3] and [4]). First of all, the in- 
spection of the phase of the wave functions of mp = 0, — 1 
components proves the formation of quantized vortices in 
these components. For mp = state the phase of the or- 
der parameter winds up by 27r (upper middle frame in 
Fig. [3]) whereas in mp = —1 component (upper right 
frame in Fig. [3|) the phase winds up by 4tt. At the same 
time the orbital angular momentum per atom equals h 
and 2h in mp = and mp = — 1 components, respec- 
tively. The circulation induced in mp = and mp = — 1 
states is, in fact, the realization of the Einstein-de Haas 
effect in cold gases. At the same time the densities ex- 
hibit non typical patterns. The density of mp = state 
consists of even number of rings when looking from a side 
(Fig. [5J upper panel) whereas in the case! of mp = — 1 
component the number of rings is odd (Fig. 2J lower 
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FIG. 1: Population of rriF = (main frame) and rriF ~ — 1 
(inset) components as a function of magnetic field (note that 
the magnetic field is pointed towards negative z-axis and only 
its values are depicted here). We observe many magnetic 
resonances corresponding to different values of magnetic field. 
Initial number of atoms in mr = 1 state is equal to N = 
5 x I0 4 and the aspect ratio is /3 = co z /uj r — 1/4. Maximal 
transfer is reached at t ~ 0.1s for each resonance. 
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FIG. 2: Population of tuf = (black) and = -1 (gray) 
components as a function of time for the first resonance. Mag- 
netic field is equal to B — — 0.12mG. 



panel). The explanation of this property will be given in 
Sec. IIV| here we just present some arguments showing 
that the number of rings in components is strictly related 
to the channel the atoms choose when they collide. 

Since the number of rings in mp = component is even 
(Fig. 21 upper panel), in the simplest case its wave func- 
tion fulfills the condition ipo( r ) ~ Yi,m, where the differ- 
ence I — vex is odd. The dipolc-dipolc interactions are small 
and can be treated as a perturbation. Therefore, the 
wave function of mp = —1 component is determined by 
nonzero amplitudes corresponding to processes in which 
only one atom changes its spin projection upon the colli- 
sion (such amplitudes are given by (Y/' im ' |^2i I^Zm)) and 
processes when both atoms change their spin projections 
(here amplitudes arc (l^',m'|^2|^m))- In the first case 
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FIG. 3: Phase (upper panel) and density (lower panel) of 
rriF = +1, 0, —1 components (left, middle, and right columns, 
respectively) in xy plane for N+i = 5 x 10 4 and j3 = 1/4. 
The presented data correspond to the first maximum in the 
resonance structure (see Fig. fTJ. Magnetic field is equal to 
B = -0.12mG. iV = 15000 and N-i = 6000 atoms both at 
their maximal transfers. 
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FIG. 4: Density of m,F = (upper panel) and mp = — 1 
(lower panel) components in xz plane for iV+i = 5 x 10 4 
atoms and f3 = 1/4. Successive frames show density patterns 
characteristic for resonances depicted in Fig. fTJ Note the 
increasing number of rings with the order of the resonance as 
well as the presence of even (odd) number of rings in tuf = 
(rriF = — 1) component. 



the amplitudes differ from zero only when m' = 1 + m 
and I' = I — 2, /, / + 2 which means that I' — m! is always 
even. Therefore, the density of = —1 state exhibits 
odd number of rings. However, in the second case V — m 
is odd what implies even number of rings in mp = — 1 
component which is not the case (see Fig. |H lower panel) . 
Evidently, the second channel is somehow closed. 
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As it was already mentioned in Section U the dipolar 
resonances can be tuned not only by the magnetic field, 
but also by the trap geometry. In order to demonstrate 
that, we have performed calculations for the spherically 
symmetric trap. Two resonances found in this case are 
shown in Fig. The properties of the first resonance 
(i.e., the one for lower magnetic field) resemble those al- 
ready discussed for axially symmetric trap. Looking, for 
example, at the densities in the xz planes we see two and 
three rings in the mp = and nip = — 1 components, 
respectively. However, the situation changes for the sec- 
ond resonance. Here, we observe inner rings structures 
for both components (see Fig. [6]). This is a signature 
that also radial excitations not only axial ones start to 
play a role. The explanation of this is given in Sec. IIVI 
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FIG. 5: Population of m,F = (main frame) and rriF ~ —1 
(inset) components as a function of magnetic field. We ob- 
serve magnetic resonances corresponding to different values 
of magnetic field. Initial number of atoms in mp = 1 state 
is equal to N+i = 5 x 10 4 and the aspect ratio is j3 = 1. 
Maximal transfer is reached at t ~ 0.3s for each resonance. 



IV. ORIGIN OF DIPOLAR RESONANCES 

To understand the origin of dipolar resonances let us 
first consider collision of two atoms, both initially in F = 
1 and mp — +1 internal state and the ground state of 
the axially symmetric harmonic oscillator, 



-(x 2 +y 2 +/3z 2 )/2 



(8) 



with corresponding single particle energy Eq = 1 + P/2. 
Here j3 = uj z /uj± is the aspect ratio of the harmonic trap. 
The oscillatory units based on the radial frequency 
are used throughout. We use a basis of single-particle 
states of harmonic oscillator and assume that contact re- 
pulsive interactions do not modify these states substan- 
tially. This assumption can be justified if occupation of 
single particle states is small. Direct evaluation of the 
two body dipolar matrix elements leads to the following 
possible outcomes of the collision: 
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FIG. 6: Phase (left column) and density (middle and right 
columns for xy and xz planes, respectively) in tuf = (upper 
panel) and tuf — — 1 (lower panel) components correspond- 
ing to the second resonance at maximal population. Initial 
number of atoms is N+i = 5 x 10 4 . Other parameters are: 
B = -0.16mG, and /3 = 1. 



i) No atom changes its spin projection. In such a case 
both atoms stay in the initial state ([5]). The elastic dipo- 
lar collisions contribute effectively to the contact inter- 
actions. 

ii) Both atoms flip their spin projection by one. To 
assure conservation of total angular momentum, the two 
atom system must gain two quanta of orbital angular 
momentum. This angular momentum can be shared by 
both atoms, then each one go to the state with mp = 
and orbital angular momentum mp = 1, or only one 
atom gains the angular momentum, i.e. it goes to the 
state with mp = and mp = 2. Second atom remains 
in mp = state. 

Evidently, in our numerical calculations we do not ob- 
serve transitions to the state with angular momentum 
mp = 2. We believe that this is due to bosonic nature of 
atoms which prefers transitions with two atoms going to 
the same final state. 

Hi) Only one atom flips its spin projection. Conserva- 
tion of total angular momentum requires that one of the 
interacting atom has to gain one quanta of angular mo- 
mentum, i.e. its wave function has mp = 1. Due to the 
symmetry of dipolar interactions the final wave function 
in z coordinate of relative motion must be antisymmetric. 

For resonances observed in our GP simulations addi- 
tional rotational quanta is shared by the atom which en- 
ters the mp = state after the collision. This can be 
checked by comparing the phases of the mp = 1 and 
mp = wave functions. There is no vortex in nif = 1 
component but there is a single vortex in mp = one, 
see Fig. [3] for example. 

To gain better understanding of observed density pro- 
files let us introduce a simplistic model where all accessi- 
ble single-particle states are limited to the lowest energy 
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FIG. 7: Spatial densities of the three lowest states coupled 
by dipolar interactions. A) ground state of the trap with 
iuf = +1; B) lowest excited state with one quanta of angular 
momentum and symmetric in the z direction. It is accessible 
when two interacting atoms flip their spins simultaneously; C) 
lowest excited state with one quanta of angular momentum 
and anti-symmetric in the z direction. It is accessible when 
one ore two interacting atoms flip their spins. 



states of the symmetry consistent with numerical results, 
i.e. to the following states 



Mr) ~ e-<* +" + ? z >/ 2 



(9a) 

<pi(r) ~ (x + iy)e- {x2+y2+Pz2)/2 (9b) 
ip 2 (r) ~ z(x + ^y)e- (a;2+ • y2+ ' 3^2)/2 (9c) 

3D graphs showing characteristic density distributions 
of all these states are shown in Fig. [7] Note that \<p 2 \ 2 
has two density rings located symmetrically with respect 
to xy plane while |<pi| 2 has only one ring in the z — 
plane. 

Treating dipolar interactions in the first order of the 
perturbation we can expand the spinor field operator in 
the basis of these single particle trap states. In a case of 
the elongated cigar shape trap the lowest excited states 
are <fi{p, 0, z) and ip 2 (fi, 0, z) and we can assume the fol- 
lowing form of spinor components of the field operator: 




<A)(?")ao 

(pi(r)ai + (p 2 (r)a 2 




(10) 



The second quantization Hamiltonian, limited by the 
expansion Eq. (|10|) consists of three parts: the single- 
particle the contact and the dipolar contributions. They 
can be written as 

H = H SP + H c + H D . (11a) 

where the single particle component is: 

H SP = (E Q + B) ftjoo + Ex a\hi + E 2 a\a 2 (lib) 

and the term proportional to B corresponds to the linear 
Zeeman shift. The contact two-body term is of the form: 

H° = Uqq ajdja a + Uu a\a\aiai + U 22 a\a\a 2 a 2 



U m aQa\a ai + Uq 2 a,Qa 2 a a 2 + U\ 2 a\a 2 aia 2 



(11c) 



and finally the two-body dipolar interaction is: 

H D = D QO alala a + ZJoiaJalaoai + D 02 alala a 2 
+ (a\a\aoao + h.c.) + T 2 (ajalao&o + h.c.) 
+ T¥\a\a\a ao + h.c.) (lid) 

All /7's, D's, and T's depend on the trap geometry and 

can be easily evaluated. Exact formulas are given in the 

( p) 

Appendix C. Parameters are dipolar integrals for 

processes when P = 1 or 2 particles are transferred from 
the initial state to the state of wave function tpt (i = 1, 2). 

Since total number of particles N = ajeto + a\ai + a\a 2 
is conserved, the Hamiltonian separates into diagonal 
blocks corresponding to a given number of particles. The 
"initial" state, 



0, 0) 



1 



(4)"lo> 



(12) 



couples, in the first order perturbation theory, with 



\N -2,2,0) 



1 



(al) N - 2 (a\f\0) 



\fW=2~j\ 



|7V-2,0,2) = 



1 



y/2(N-2)\ 



(13a) 
(13b) 
(13c) 



where |0) is the vacuum. Let us notice that dipolar ma- 
trix element corresponding to two particles going to tp 2 (r) 

state is proportional to ^/2N{N — 1) T^p which in the 

limit of large N scales as ~ NT 2 ■ The corresponding 
matrix element associated with only one particle jumping 
to this state is \fN (N—l)^ 1 ^ and scales with number of 
particles as ~ N^^T^" ■ Evidently in the limit of large N 
atoms enter the state ^(r) one by one - only one atom 
per every two-body dipolar collision flips its spin. The 

processes driven by term proportional to dominate 
the collisions driven by the for large N . 

For N — 2 the four states listed above span the whole 
subspace. The spectrum obtained from the exact diag- 
onalization is presented in Fig. [8l The figure explains 
origin of resonant behavior of the system. For B = 
energy of each state is a sum of kinetic energy, poten- 
tial energy, and contact interaction energy. Energy re- 
lated to the dipole-dipole interaction is very small and 
can be ignored. When the magnetic field is switched 
on the additional Zeeman energy must be taken into ac- 
count for atoms which are in states with mf = +1. As 
we see for some magnetic fields the energies of different 
states become equal and then dipolar interactions start 
to play a crucial role. First, they cancel degeneracy of 
the states, second, they can lead to noticeable transfer 
between them. This is how the dipolar resonances do 
appear. 
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Three avoided crossings seen in the Fig. [8] correspond 
to three different resonances. The width of every res- 
onance is determined by the magnitude of the dipolar 
energy and is very small. For two rubidium atoms in a 
magnetic trap of frequency ~ 100Hz it is of the order of 
AB ~ 10~ 4 mG. Such a precision in control of the mag- 
netic field is extremely difficult - beyond experimental 
reach. However, the dipolar energy grows with number 
of atoms as discussed above. Therefore for N ~ 10 4 (as 
studied in section IIII[) , widths of resonances increase by 
several orders of magnitude AB ~ 10 _1 mG giving a hope 
for an experimental realization. Thus for large systems 
the effects discussed could be observed in a rubidium con- 
densate in experiments performed at ultra-low magnetic 
fields requiring a special shielding of all external mag- 
netic fields. On the other hand, large atomic densities 
can be reached in optical lattices with few atoms per site 
only. In such a situation chromium atoms are preferred 
because of relatively large magnetic moment. For typical 
optical lattices the width of resonances for two chromium 
atoms in a single site is about AB ~ 10~ 1 mG. Moreover, 
in the case of optical latice, the width of the resonance 
is significantly increased due to the finite energy width 
of excited bands as seen in experiment pol ]. Again spe- 
cial shielding is required for experimental observation of 
resonant dipolar interactions. If the magnetic field docs 
not match a resonant value, the dipolar interactions can 
be totally ignored. 



Weakness of dipolar interactions is responsible for a 
very small width of the resonances making them difficult 
to detect. On the other hand weak dipolar interactions 
have one very important advantage: all the resonances 
shown in Fig. [8] are well separated. Choosing an appro- 
priate value of the magnetic field a particular resonance 
can be addressed, and dipolar interactions can be tuned 
to couple desired center of mass states. Obviously far 
from all the resonances the dipolar interaction can be set 
to zero. 
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FIG. 8: Energy of two-atom quantum states as a function of 
external magnetic field for ui z /u)±_ — 1/4. For some values of 
external magnetic field energies of some states become equal 
due to the linear Zeeman effect. Two insets show real and 
avoiding crossings. Note that avoiding crossings appear only 
for states coupled by dipolar interactions. 



V. CONCLUSIONS 

In conclusion, we have shown the existence of dipo- 
lar resonances in rubidium spinor condensates. The 
resonances occur when the Zeeman energy of atoms in 
m,F = 1 component matches the excitation energy of 
a spatial state which is populated upon the transfer. 
Symmetries of dipolar interactions force the atoms in 
mp = 0,-1 states to circulate around the quantization 
axis. In fact, singly and doubly quantized vortices are 
formed in mp = 0,-1 components, respectively. Since 
dipolar interactions are weak the resonances are narrow. 
It means that dipolar resonances are very selective, i.e. 
one can choose the transition to the proper state by 
choosing the magnetic field. Hence, dipolar resonances 
seem to be a route to the observation of the Einstein-de 
Haas effect, as well as other phenomena related to dipolar 
interactions, in weak dipolar systems. 



Exact value of the resonant magnetic field depends on 
energy difference of resonantly coupled states at B = 
and it is a linear function of number of atoms in the initial 
state. For realistic values of scattering lengths the con- 
tact interaction in the initial \N, 0, 0) state grows with N 
faster then the contact interaction in three other consid- 
ered here states. Therefore, the total energy of this state 
can be larger then the energy of the resonantly coupled 
state. It clearly follows from Fig. [8] that for large N the 
corresponding resonant magnetic field is negative, i.e. its 
direction is opposite to the direction of the magnetic mo- 
ment of atoms. This simple observation explains why in 
our numerical calculations discussed in sec. Mil a direc- 
tion of the magnetic field was inverted with respect to the 
direction used for the preparation of the ground state. 
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Appendix A: Equation of motion 



and 



The equation of motion ([3]) is obtained from the 
Heisenberg equations for the field operators ipi{r) for the 
system described by the Hamiltonian, H (Q}. Since the 
Hamiltonian ([T]) splits in single-particle part and two- 
particle parts related to contact and dipolar interactions, 
the commutators [4>i(r), H] as well as the right hand-site 
of Eq. §5§ get the form of the sum of three terms of differ- 
ent physical origin. The single-particle part is diagonal 
and given by H sp = -^V 2 + V ext - fim F B. 

The term T-L c in Eq. JSj, which is related to the contact 
interactions has the following diagonal elements 

%li = (c + c 2 ) Tp\ipi + (c + c 2 ) Vwo 

+ (co- c 2 )V>l 1 V'_i 
^00 = ( c + c 2) + c %ipo 

+ (co + c 2 ) V'LiV'-i 

= (co - c 2 ) ^I^i + (cq + C 2 ) ^J^Q 



+ (c + c 2 )^L 1 ^_j 



(Al) 



These elements describe the collisions of atoms that pre- 
serve the projection of the spin for each atom. The off- 
diagonal elements, on the other hand, are responsible for 
collisions changing separately the atomic spin projections 
but conserving the projection of total spin. They yield 









^0-1 


= c 2 %ipi 




T-Ci-i 


= 0. 


(A2) 



For dipolar interactions, H d , in Eq. ([3]) one has 



nUr) 



d 3 r'$t(r')Vf jM (r-r%(r'), (A3) 



where 



V l i u {r-r') = 



\r — r 



'13 



FijFki 



-[F^r-r'^Fuir-r')}. (A4) 



Only two elements oiH d matrix are independent, namely 



i 



,(*-*') 



l\2 



r — r 



\r — r 



/lo 



-3^= / d 3 -' 
\/2 



x(Wo + W-i) 



j-^—^KZ-zo-^-J,') 



'71 

fit,?, 



3„/ 



r — r 



— [(a- -a/) +*(»-»') 



(A5) 



H d w (r) = -3^jd 



3 r / [(s-s')-i(2/-y')](z-z' 

I — i**' 1 5 



"2^ 



d 3 r 



, [(s-s / )-*(y-y , )] 2 ,7t r . r 1 7 v 

i 777 (W0+ IpO 1p-l) 



1 



r — r 



/|3 



For other elements one has 



r — r' 

3 (x - x'f + (y - iff 
2 |r-r'| 5 



H 



o-i 



'MO) n -l-l 



njd 



nid _ nia 
H 1 _ 1 — H 



(A6) 



(A7) 



The Hfj terms (both diagonal and off-diagonal) are 
responsible for the change of the total spin projection of 
colliding atoms. For example, even if components mp = 
and mp = — 1 are initially empty then due to the first 
part in H d there exists a process when one of two atoms 
being initially in mp = 1 state populates mp = state 
after collision. 



Appendix B: Numerical details 

To solve Eq. ([3]) we neglect the quantum fluctuations 
and replace the field operator V>j(r) (V>f (t*)) by an order 
parameter ipi( r ) (V'i*( r )) for each component. Next, we 
split the right-hand side operator of (J3)) in the following 
way 



n = w 



U d = T+V 



(Bl) 



where T is the kinetic energy operator which is a diagonal 
matrix and V represents the potential energy including 
the trapping potential, the interaction with the external 
magnetic field, and the contact and dipolar interactions. 
Now, we apply the split-operator method to solve the 
mean- field version of Eq. (J3)). First, we calculate the 
evolution of an order parameter according to the poten- 
tial energy. To this end, we diagonalize the matrix V 
at each spatial point V = PDP^ 1 , where the matrix P 
consists of the eigenvectors written in columns and the 
matrix D is diagonal with the eigenvalues on the diago- 
nal. The evolution by the time step At is given by 




= p e - i£)A */fip-i 




(B2) 



and is calculated in three steps: first, the spinor or- 
der parameter is multiplied by the matrix next by 
e -iDAt/h^ wn j cn j s a diagonal matrix since D is diagonal, 
and finally by the matrix P. Then we calculate the evo- 
lution due to the kinetic part. Since the kinetic energy 
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operator T is diagonal wc apply the Fourier technique to 
each component in a way as it is usually done in the case 
of a scalar Gross-Pitacvskii equation. 

All integrals appearing in formulas (|A5[) and (| A6|l are 
the convolutions and we use the Fourier transform tech- 
nique to calculate them. The Fourier transform of the 
convolution of two functions is the product of the Fourier 
transforms of these functions. Therefore, it is useful to 
have analytical formulas for the Fourier transforms of the 
components of the convolutions that do not change dur- 
ing the condensate evolution. The appropriate Fourier 
transforms are found in the following way. First, two ro- 
tations of coordinates system (by angles: /3 around z-axis 
and a around y-axis) are applied to make the vector k - 
the argument of the Fourier transform - parallel to the z- 
axis and next, the regularization procedure as described 
in Ref. 1241 is used. The result reads 



T 
T 
T 
T 
T 



1 



- 3, 



= —-^-(1 — 3 cos 2 a) 



1 3 x" + y 
rP ~ 2~ \r 



2 1 



2?r ft o 2 \ 

= -^-(1 — 3 cos a) 



(x - iy)z 



\r\° 
(x + iy)z 



(x - iyf 



= ^L e -^ sin 2a 
3 

= — e lP sin 2a 
3 

= -^ e -^sin 2 a 



(B3) 



with the angles a and /3 defined as follows: cos a = k z /k 
and sin j3 = k y j ^Jk 2 — k 2 . 

Even with a help of formulas (|B3j) the evolution of 
mean-field version of Eq. ([3]) is a time-consuming com- 
putational task. We would like to point out that each 
time step consists of taking nine numerical Fast Fourier 
Transforms (three for determining the evolution of the 
spinor wave function governed by the kinetic energy and 
six for calculating T~L d matrix) and also doing the diago- 
nalization of matrix V at each spatial point. Therefore 
the parallel computing with OpenMP directives has been 
used in our calculations. 

Size of the grid was adjusted to the trap geometry 
- we have chosen 2 5 x 2 5 x 2 5 grid points in the case 
of spherically symmetric traps (with the spatial steps 
equal to Ax = Ay = Az = 0.55 osc. units, where osc. 
miit= y^h/muj) and 2 5 x 2 5 x 2 6 grid points for cigar-like 
traps (with the spatial steps equal to Ax = Ay = 0.275 
and Az = 0.375 osc. units). The time step was chosen to 
At = 2.5 x 10~ 5 osc. units, where osc. unit=l/w. Part 
of the results was recomputed with twice smaller spatial 
step (two times bigger grid size in every dimension) as a 
self-test. 

The ground state was obtained by evolving the mean- 
field version of Eq. Q in imaginary time at the value 
of magnetic field B = 0.73 mG. Such magnetic field was 
large enough to result in having almost all atoms inmj? = 
1 Zeeman sublevel. 

Appendix C: Two-body matrix elements 



defined by ©. 
Contact integrals: 



£/nn 


co + 
2 


c 2 / 
J 


d 3 r \ifin(r)\ 4 




(Cla) 


u n 


co j 
~~ 2 J 


d 3 r 1 


g3i(r)| 4 




(Clb) 


T J 


CO i 

= Ti 


Q T \ 






{HI r-\ 

(Clc) 


^01 






/d 3 r \wn(r)\ 2 
f u ' I ^0 \ 1 ) 1 


09i <r)\ 2 




U02 


= (co 4 


-c 2 ) 


[d 3 r Mr)| 2 




(Cle) 


U12 




fd 3 r 


ipi(r)\ 2 \ip 2 (r 


I 2 


(Clf) 




co t 
" 2 J 


d 3 r 


1 2 




(Clg) 



Dipolar interactions which do not change the spin of 
atoms: 



Ao = 



^ / /d 3 rdV \Mr)\ 2 .h(r-r')\Mr')\ 2 



(C2a) 



Ai = - 



A2 = - 



d 3 rdV v>;(r)vJo(r)/i(r - r') V l(r')<p Q {r') 

(C2b) 

d 3 rdV ^(r)v?o(r)/i(r - r>|(r')w>(r') 

(C2c) 



Dipolar interactions which change the spin of the 
atoms: 

(2) _ V* I 1 , 3 „ j3„/ */ 



Tl = 2 J J d rd r ' <P*i( r )Mr)f2(r- r')^{r')^{r') 

(C3a) 

T 2 (1) =M 2 / /d 3 rdV <pl{r)y o {r)h{r-r')<pl{r')y {r') 

(C3b) 

if 5 = yJ /d^dV ^(r)^ (r)/2(r- r>5(r')po(r') 

(C3c) 



Functions fi have a form: 

l-'Sn 2 



fi{r) 
Mr) 
Mr) 



3 (n x - in y ) 2 
2 ^ 

3V2 n z (n x + in y ) 



(C4a) 
(C4b) 
(C4c) 



In this section we present explicit formulas for two- 
body matrix elements of the Hamiltonian ([1]) in the basis 



10 



A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005). [15 
Q. Beaufils, R. Chicireami, T. Zanon, B. Laburthe-Tolra, 
E. Marechal, L. Vernac, J.-C. Keller, and O. Gorceix, [16 
Phys. Rev. A 77, 061601 (2008). 

T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. [17 
Pfau, Rep. Prog. Phys. 72, 1 26401 (2009). 

M. Ueda and Y. Kawaguchi, larXiv: 1001.20721 (2010) . [18 
T. Lahaye, T. Koch, B. Frohlich, M. Fattori, J. Metz, A. 
Griesmaier, S. Giovanazzi, and T. Pfau, Nature 448, 672 [19 
(2007). 

T. Lahaye, J. Metz, B. Frohlich, T. Koch, M. Meister, 

A. Griesmaier, T. Pfau, H. Saito, Y. Kawaguchi, and M. [20 
Ueda, Phys. Rev. Lett. 101, 080401 (2008). 
S. Yi and H. Pu, Phys. Rev. Lett. 97, 020401 (2006). 
Y. Kawaguchi, H. Saito, and M. Ueda, Phys. Rev. Lett. [21 
98, 110406 (2007). 

K. Gawryluk, M. Brewczyk, K. Bongs, and M. Gajda, [22 
Phys. Rev. Lett. 99, 130401 (2007). 

B. Sun and L. You, Phys. Rev. Lett. 99, 150402 (2007). 
J.-N. Zhang, L. He, H. Pu, C.-P. Sun, and S. Yi, Phys. [23 
Rev. A 79, 033615 (2009). 

S. Hoshi and H. Saito, Phys. Rev. A 81, 013627 (2010). 
T. Swislocki, M. Brewczyk, M. Gajda, and K. Rzazewski, [24 
Phys. Rev. A 81, 033604 (2010). [25 
[14] M. Vengalattore, S.R. Leslie, J. Guzman, and D.M. 



Stamper-Kurn, Phys. Rev. Lett. 100, 170403 (2008). 
Y. Kawaguchi, H. Saito, K. Kudo, and M. Ueda, Phys. 
Rev. A 82, 043627 (2010). 

A. Einstein and W.J. de Haas, Verh. Dtsch. Phys. Ges. 
17, 152 (1915). 

Y. Kawaguchi, H. Saito, and M. Ueda, Phys. Rev. Lett 
96, 080405 (2006). 

L. Santos and T. Pfau, Phys. Rev. Lett. 96, 190404 
(2006). 

B. Pasquiou, G. Bismut, A. Crubellier, E. Marechal, P. 
Pedri, L. Vernac, O. Gorceix, and B. Laburthe-Tolra, 
Phys. Rev. A 81, 042716 (2010). 

B. Pasquiou, G. Bismut, E. Marechal, P. Pedri, L. 
Vernac, O. Gorceix, and B. Laburthe-Tolra, Phys. Rev. 
Lett. 106, 015301 (2011). 

T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998); T. Ohmi and 
K. Machida, J. Phys. Soc. Jap. 67, 1822 (1998). 
E.G.M. van Kempen, S.J.J.M.F. Kokkelmans, D.J. 
Heinzen, and B.J. Verhaar, Phys. Rev. Lett. 88, 093201 
(2002). 

C. J. Pethick and H. Smith, Bose-Einstein Condensa- 
tion in Dilute Gases (Cambridge University Press, Cam- 
bridge, United Kingdom, 2002). 

K. Goral and L. Santos, Phys. Rev. A 66, 023613 (2002). 
B. Pasquiou, et. al., Phys. Rev. Lett. 106, 015301 (2011). 



